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ABSTRACT 

We use early-time photometry and spectroscopy of 12 Type II plateau supernovae 
(SNe IIP) to derive their distances using the expanding photosphere method (EPM). 
We perform this study using two sets of Type II supernova (SN II) atmosphere models, 
three filter subsets {{BV}, {BVI}, {VI}), and two methods for the host-galaxy ex- 
tinction, which leads to 12 Hubble diagrams. We find that systematic differences in the 
atmosphere models lead to ~ 50% differences in the EPM distances and to a value of Hq 
between 52 and 101 km s~^ Mpc~^. Using the {VI} filter subset we obtain the lowest 
dispersion in the Hubble diagram, cj^ = 0.32 mag. We also apply the EPM analysis 
to the well-observed SN IIP 1999em. With the {VI} filter subset we derive a distance 
ranging from 9.3 it 0.5 Mpc to 13.9 it 1.4 Mpc depending on the atmosphere model 
employed. 

Subject headings: (stars:) supernovae: general - galaxies: distances and redshifts 

1. Introduction 

Type II supernovae (SNe II) are understood as the result of the final gravitational collapse 
of massive stars {M > 8 Mq) that at the moment of the explosion have most of their hydrogen 
envelope intact. The energy released in the explosion is typically ~ 10^^ erg (mainly radiated in 
the form of neutrinos), and the luminosity of the SN during the first few months after explosion can 
be comparable to the total luminosity of its host galaxy. These objects have been classified based 
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on their light curves into Type IIP (plateau) and Type IIL (linear) (e.g., iPatat et al.lll994l ). The 
former present a nearly constant optical luminosity during the photospheric phase (~ 100 days after 
explosion), while the latter show a slow decline in luminosity during that phase. However, there are 
some SN II events, such as the SN 1987A, that show peculiar photometric properties. Also, studies 
of SN II spectra have reveale d the existence of a subclass, char acterized by the presence of narrow 

1991al lbl^ 



spectral lines, called SNe Iln (jSchlegel 



1990; 



Filippenko 



which are mo st likely originated 
from the interaction of the SN ejecta with pre-existing circumstellar material; see iFilippenkd (jl997l ) 
for a general review of SN spectra. 

Due to their high intrinsic luminosities, SNe II have great potential as extragalactic distance 
indicators. To date, several methods have been proposed to derive dist ances to SNe II, but two 
are the most commonly used: the expanding photosphere method (E PM) (jKirshner &: Kwanlll974l ) 
and the standardized candle method (SCM) (iHamuv &: Pintd |2002| ) . The former is a geometrical 
technique that relates the photospheric radius and the angular radius of a SN in ord er to derive its 



distance, and has been applied to several SNe to derive the Hubble constant (e.g., ISchmidt et al 



19921 ). The EPM is independent of the extragalactic distance ladder, and therefore does not need 



any external calibration. The SCM is based on the observed relation between expansion velocity 
and l uminosity of SNe I IP. Recently, this method has been applied to a sample of high-redshift 
SNe (iNugent et al.l 120061 ). Other methods have also been used to det ermine dista,nces t o SNe II, 
such as the spectral-fitting expa nding atmo s phere method (SEAM) (| Baron et al.l |200J) and the 
plateau-tail relation proposed by iNadvozhinI (|2003l ) . 

In this work we apply the EPM using early spectroscopy and photometry of 12 SNe IIP in 



order to derive their distances. We ap 



ply the method using two sets of SN II atmosphere models 



([Eastman et al.lll996l : lDessart &: Hillier 



2005a|), three filter subsets {{BV}, {BVI}, {VI}), and two 



methods for the host-galaxy extinction. The different combinations lead to 12 Hubble diagrams. 
Section 2 of this paper describes the photometric and spectroscopic observations. In § 3, the EPM 
is presented, and we apply it to 12 SNe IIP. The results are discussed in § 4. We compare our EPM 
distances with results from other methods and with previous EPM analyses. We also discuss the 
error analysis and the effect of reddening on the EPM distances. We show 12 Hubble diagrams and 
the corresponding Hubble constants, and we propose an external calibration for the EPM. Finally, 
we summarize our conclusions in S 5. 



Observations 



In this work we use photometry and spectroscopy from four SN follow-up programs: the Cerro 
Tololo supernova program (1986-1996), the Calan/Tololo supernova survey (CT; 1990-1993), the 
Supernova Optical and Infrared Survey (SOIRS; 1999-2000) and the Carnegie Type II Supernova 
Program (CATS; 2002-2003). During these programs optical (and some IR) photometry and spec- 
troscopy were obtained for nearly 100 SNe, 51 of which belong to the Type II class. All of the 
optical data have already been reduced and will soon be published (iHamuy et al.l 120081 ). We also 
complemented our dataset with some spectra from various coauthors of this paper. 
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2.1. Photometry 



Direct images of SNe were obtained with telescopes from four different observatories: the Cerro 
Tololo Inter- American Observatory (CTIO), the Las Campanas Observatory (LCO), the European 
Southern Observatory (ESO) in La Siha, and the Steward Observatory (SO). Several telescopes 
and instruments were used to obtain the photometry, which is listed in an el ectronic Table. In all 



cases CCD de tectors and standard Johnson-Kron-Cousins UB Fi? 72^ filters (| Johnson et al.l 119661 : 



Cousinslll97lh were employed. For a small subset of SNe observations in the JiJiT filters were also 



ob tained. The data re duction was performed using according to the procedure described 

by lHamuy et al.l (120081 ). The error in the photometry ranges between 0.01-0.06 mag, with a typical 
value of 0.02 mag. 

The optical light curves of all the SNe used in this work are shown in Figures [iHSl clearly 
revealing the plateau nature of all these events. 



2.2. Spectroscopy 



Low resolution {R ~ 1000) optical spectra (wavelength range ~ 3200-10000 A) were taken for 
each SN at various epochs using telescopes and instruments from four different observatories. An 
electronic Table lists all the telescopes and instruments used for the spectroscopy. Most of the 
spectra were obtained with the slit along the parallactic angle (Filippenko 1982). The wavelength 
calibration was performed using comparison-lamp spectra taken at the position of each SN . The 



flux calibration was done via observations of flux-standard stars (Hamuv et al. 



1992 



199J). For 



more details on the observational procedures see iHamuv et al.l (|2008l ). 

The spectra were shifted to the rest frame using the heliocentric redshifts given in Table [1] in 
order to measure the SN ejecta velocities. In seven cases we were able to measure the redshifts 
from narrow emission lines of H II regio ns at the SN positio n (see Table [1]). Also, in one case 
(SN 1999em) we adopted the value from [Leonard et al.l (2002) which corresponds to the redshift 
measured at the SN position. In four cases we were unable to extract this information from our 
data, and we had to rely on redshifts of the host-galaxy nuclei; this does not take into account the 
rotation velocities of the host galaxies, which are typically v ~ 200 km s~^. 



2.3. Sample of Supernovae Used in this Work 

Fifty-one SNe II were observed in the surveys described above. We cut this sample according 
to the EPM requirements, which are (1) the optical SN light curve (V and / bands) must show a 
nearly constant luminosity during the photospheric phase, i.e, the SN must belong to the SN IIP 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
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class (see Figures [TH3|) ; (2) the SN must have early-time photometry; and (3) the SN must have at 
least three early spectroscopic observations. The necessity for all of these requirements, discussed 
in § 3.6, reduced the sample to only 11 SNe. We also adde d the SN IIP 1999gi to our sample, which 
has extensive photometry and spectroscopy published by [Leonard et al.l (j2002bl ). 



3. The Expanding Photosphere Method 

3.1. Basic Principles 

The EPM is a geometrical technique that relates an angular size and a physical size of a SN, 
in order to derive its distance. Although the angular radius of a SN cannot be resolved spatially 
with current optical instrumentation, it can be derived assuming a spherically symm etric expanding 



phot osphere (a reasonable assumption for SNe IIP at early times, as discussed by [Leonard et al 
that radiates as a black body "diluted" by a factor <^^. Specifically, 



R 
D 



{l + z)fx 



<2,5y(r)10-0-4[^W+A'{A') 



(1) 



where R is the photospheric radius, D is the distance to the SN, fx is the observed flux density, A is 
the observed wavelength, B^i is the Planck function in the SN rest frame, T is the color temperature, 
A = A/(l + z) is the corresponding wavelength in the SN rest frame, ^(A) is the foreground dust 
extinction and A (A ) is the host-galaxy extinction. The factor ("^/ (known as "distance correction 
factor" or "dilution factor") accounts for the fact that a SN does not radiate as a perfect black 
body; there is flux dilution caused by grey electron scattering which makes the photosphere (defined 
as the region of total optical depth r = 2/3) form in a layer above the thermalization surface. Also, 
the dilution factor accounts for line blanketing in the SN atmosphere. Since electron scattering is 
the main source of continuum opacity, the total opacity is essentially grey, and the photospheri c 
angular radius is independent of wavelength in the optical and near- infrared ([Eastman et al.lll996l ). 
which explains why R and 9 do not carry a wavelength subscript. 

Because the gravitational binding energy (U ~ 10^^ erg) of a SN progenitor is far less than the 
expansion kinetic energy {E ~ 10^^ erg) of the ejecta, it is reasonable to assume free expansion. 
This assumption is supported by hydro dynamical models which show that the different layers of 
the ejecta reach ~ 95% of their terminal velocities ~ 1 day after the explosion. During this brief 
period the re is a transition from an ac celeration phase due to the SN explosion, to homologous 
expansion ( UtrobirJ [2007[ : [BerstenI [20081 ). Due to the high expansion velocities (~ 10000 km s~^), 
the initial radius (typically Rq ~ 10^^ cm for a red supergiant) can be neglected after ~ 1 day from 
explosion; hence after that period the physical radius of the SN can be approximated by 



(2) 
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where v is the photospheric velocity and to is the explosion date. Combining (1) and (2) we obtain 

^ _ (tj - to) 

where 9i and Vi are the derived quantities measured at time tj, which are estimated following the 
steps explained in the next sections. Equation 3 shows that the quantity 6/v increases linearly with 
time, so D and tg can be derived from at least two spectroscopic and photometric observations. 
More observations allow us to check the internal consistency of the method. 



3.2. Dilution Factors 

The dilution factors correspond to the ratio of the luminosity of a SN atmosphere model {L^r ) 
and the corresponding black-body luminosity. 



^2 



(4) 



In practice, the dilution factors must be derived for the same filter subsets employed to determine 
the color temperature (T) of a SN. In this work we focus on three differen t optical filter subset s 
{{BV}, {BVI} and I V^Jj), and we us e d two S N atmosphere models, those bv lEastman et al.l (jl996l ) 
(E9 6 hereafter) andlDessart fc Hillieii (j2005bl ) (D05 hereafter), to compute the dilution factors. See 
also lDessart fc HiUierl (j2005al ) for more details of the input parameters of the D05 models. Because 
the color temperature of each SN was determined from colors measured in the observer's rest frame, 
both the atmosphere models and the black-body function must be redshifted; thus, the dilution 
factors must be computed for the specific redshift of each SN. 

We computed B, F,/ synthetic magnitudes using 58 spectra from E96 atmosphere models and 
138 spectra from D05 atmosphere models. For each filter subset S (that is, S = {BV}^ {BVI}^ {VI}), 
we fit black-body functions in the SN rest frame B^'{Ts), and solved for and (^s,z by minimizing 
the quantity 



^[M^ + 51og( 
AGS 



R 



10 pc 



+ 51og(C5,.)-^(T.,z)]^ 



(5) 



Here R is the photospheric radius, Mj is the redshifted synthetic absolute magnitude of the at- 
mosphere model for a band with central wavelength A, and b-^{Ts, z) is the synthetic magnitude of 

vrSy (T,)10-0-4[^W+^'(^')l/(l + z), given by 
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h = -2.5 log.o / ^ ,,(1 + ,) S{X)dX + ZP, (6) 

where S{\) is the fi lter transmision function and ZP is the zero point of the photometric system 



(jHamuv et al.l |2001| ). The constants h and c are the Planck constant and the speed of light, 
respectively. Clearly, the dilution factors depend on the specific redshift of the SN and on the filter 
subset used to obtain temperature of the models. Figure H] shows the resulting dilution factors 
versus temperature at z = 0. We performed polynomial fits to C,{Ts) of the form 

aTs) =j:b.A^^\ . (7) 

Table [2] lists the hgj coefficients at z = for three filter subsets and both atmosphere models (E96 
and DOS). The corresponding polynomial fits are shown in Figure HI 

The DOS dilution factors are quite insensitive to the color temperature above ~ 9000 K, and 
lie around 0.5, while at lower temperatures they increase sharply with decreasing temperature, 
reaching a value over unity below ~ SOOO K. The E96 dilution factors present the same pattern, but 
they are sy stematically lower than th e DOS dilution factors by ~ IS %. The origin of these differences 



is unclear. iDessart &: Hillied (|2005al ) discuss that the discrepancy might be related to the different 
approach used to handle relativistic terms. Also, DOS solved the non-LTE (local thermodynamic 
equilibrium) (non-LTE) problem for all the species, and employed a very complex atom model. 
E96, on the other hand, solved the non-LTE problem for a few species, while for the rest of the 
metals the excitation and ionization were assumed to be given by the Saha-Boltzmann equation, 
and the opacity was taken as pure scattering. Another important difference between the E96 and 
DOS dilution factors is the dependence on the parameters involved in the atmosphere modelling. 
While the E96 dilution factors show little sensivity to a broad range of phyical parameters other 
than temperature, the DOS models show a larger dispersion at a given color temperature. However, 
this is also due to DOS models covering a larger range of radii, density profiles {p oc r^"^) and 
temperature than E96. On average, the E96 models lead to a dispersion of cj ~ 0.03 in while 
the DOS models yield a ~ 0.07. 



3.3. Angular Radii 

An apparent angular radius (^Cs) s-nd a color temperature (T^) of the SN can be obtained 
by fitting a Planck function Byi{Ts) to the observed broad-band magnitudes (see eq. 1). Here 
iS* is the filter subset combination, i.e., S = {BV}, {BVI}, {VI}. Since we have two unknowns 
{6Cs-iTs), the subsets must contain at least two filters. In order to derive these parameters, we 
used a least-squares technique at each spectroscopic observation epoch (see §3.6) by minimizing 
the quantity 



-7- 



2 ^ [mj + 5log{eCs,z)-hiTs,z)Y 

X -2^ -2 • (»j 

s 

Here, is the apparent magnitude in the filter with central wavelength A (i.e., G /}), 
dm is the photometric error in the magnitude mj, and is defined in eq. 6. Because Cs is mainly 
a function of the color temperature (Fig. S]), it is possible to use to solve for C^g and determine 
the true angular radius d from QC^g. 



3.4. Physical Radii 



Once d is determined, the next step is to measure the photospheric velocity (see eq. 3). The 
photospheric velocity of the SN at a given epoch can be obtained from the absorption lines in the 
spectra. We measured velocitie^ from the minima of Hq, H/3, H7, and Fe 11 A 5169 lines, for all 
12 SNe. The tables that list the spectroscopic velocities are available in electronic format. Figures 
[5][7]show the temporal evolution of the spectral line velocities. 

To date the photospheric velocities have been e stimated using weak s pectral absorption fe atures 
such as Fe 11 AA5169, 5018, 4924, and Sc 11 A4670 (jSchmidt et al.lll992l : Eeonard et ahlbood l. The 
physical assumption is that these lines are weak and formed near the photosphere of the SN. 
However, there are two problems with this approach: (1) at early times the spectra are dominated 
by Balmer lines and the weak lines are absent, and (2) the sy nthetic spectra show tha t even the 
weak lines do not necessarily yield true photospheric velocities (jPessart Hillierl l2006l ) . One way 
to circumvent these problems is to use the Balmer lines which are present in the spectra over most 
of the evolution o f the SN. Although the Balmer lines are optically thicker than the Fe 11 lines. 



Dessart &: Hillierl (120061 ) argue that, contrary to what is usually believed, optically thick lines do 



not necessarily overestimate the photospheric velocity, and the offset from the photospheric velocity 
can be measured from the synthetic spectra. In this work we decided to use the minimum of the 
H/3 absorption line to derive the photospheric velocity because this line is present during the entire 
plateau phase, it can be easily identified, and it does not present any blend, at least in the first 
~ 50 days after explosion. 

To convert from observed H/3 spectroscopic velocities to true photospheric velocities we used 
the synthetic spectra from E96 and DOS. Figure [8] shows (in red) the ratio of H/3 velocity and the 
photospheric velocity, as a function of H/3 velocity for all of the D05 models. Note that the DOS 
models predict that the H/3 line forms quite close to the photosphere at all epochs (for all values 
of fH/3). Also plotted in Figure [8] (in blue) are the E96 models which confirm that the H/3 forms 



^We employed the non-relativistic formulae to derive the expansion velocities from the doppler shift of the ab- 
sorption lines. This approach is reasonable because the highest velocity used in this work is ~0.035c, for which 
the difference in the velocity using the relativistic and non-relativistic formulae is less than 2%. However, typical 
velocities employed in the EPM analysis are of 0.02c, for which the difference is ~ 1%. 
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close to the photosphere at early epochs, when vnp is high. However, at later epochs (lower vup) 
E96 predict that H/3 forms in outer layers (higher velocities) than D05. It is also important to note 
that the E96 models cover a shorter range in velocity (~ 4500-12000 km s"^) than the D05 models 
(~ 2000-17000 km s~^), which restricts the EPM analysis using the E96 models. 

To derive the ratio between the H/3 and the photospheric velocity we used a polynomial fit of 
the form 

^ = jz^AvHpy^ (9) 

Uphot 

(see Figure [8|). The aj coefficients are listed in Table [3j The E96 models lead to a dispersion of 
a = 0.06 and the DOS models to cr = 0.04. The photospheric velocity Vi can be obtained from a 
measurement of vup'- 

- - . (10) 



2 

j=0 



In order to examine which of the adopted photospheric velocity conversions was closer to 
reality, we compared the ratio between the Ha and H/3 velocities measured from the observed 
spectra of our sample of SNe and from the synthetic spectra of the E96 and DOS models. Figure 
[9] shows the Yia/Hj3 velocity ratio as a function of the H/? velocity. It can be seen that, while 
there is good agreement between theory and observations at high H/3 velocities (~ 7000-10500 
km s~^), the DOS models underestimate the Ha velocities (or overestimate the H/? velocities) at 
lower expansion velocities. This could be due to time-dependeri ce effects that are not include in 



the DOS models that assume steady-state (jPessart fc Hillierll2008l ). On the other hand, the Ha/H/3 



velocity ratio predicted by the E96 models is in good agreement with the observations at all H/3 
velocities, although there are few models below ~ 6000 km s~^ to draw strong conclusions. This 
suggests that E96 predict more realistic line profiles in the SN ejecta than DOS and therefore should 
provide a better photospheric velocity conversion. 



3.5. Extinction 



To estimate the amount of Galactic foreground extinction we used the IR dust maps of 
Schlegel et al.l (|l998l ). Table H] summarizes the foreground extinction adopted. In this work we 
used two different methods for the determination of host-galaxy reddenings of our SN sample: a 
spectroscopic method (DES hereafter), and a method based on the color evolution of the SNe 
(OLI hereafter). The former consists in fitting different model spectra to the early-time spectra 
of a S N. The two fitting parameters are the amo unt of reddening and the photospheric temper- 



ature (Dessart &: Hillier 



20061 : iDessart et al.l l2008l ) . The color-based technique was developed by 
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Olivares et alj (|2008l ) and is based on the assumption th at the color at the end of the plateau phase 



is the same for all SNe IIP. In both cases we adopted the lCardelli et al.1 (|l989l ) extinction law (with 
Rv = 3.1). 

Table U] lists the host galaxy visual extinction values Ay obtained from both methods. Also, 
Figure [10] shows the OLI versus DES visual extinctions. As can be seen, there are no systematic 
differences between the models. However, there are individual differences, especially in five SNe, 
whose names are explicitly marked in the plot. 



3.6. Implementation of EPM 

The EPM method is only valid in the optically thick phase of a H-rich expanding atmosphere. 
Observationally, this period corresponds to the plateau phase of Type II SNe and thus justifies our 
first selection criterion in § 2.3. 

The EPM requires at least two simultaneous photometric and spectroscopic observations (see 
eq. 3), but we recommend the use of at least three points in order to obtain an internal check. The 
photometry is used to determine the angular size of the SN and the spectroscopy is used to measure 
the expansion velocities of the SN. The requirement of simultaneous photometric and spectroscopic 
observations is usually not accomplished; the photometry and the spectroscopy of a SN are taken 
at different epochs. To overcome this problem, it is necessary to interpolate the photometry or the 
velocities measured from the spectra. In this work we decided to interpolate the photometry for 
two reasons: (1) the number of photometric observations in our sample of SNe is far greater than 
the number of spectroscopic observations, and (2) the optical apparent magnitude of SNe II-P is 
nearly constant during the plateau phase, making the photometry interpolation more reliable than 
the velocity interpolation, which has a steeper dependence with time. To interpolate a magnitude 
at the epoch of a given spectroscopic observation, we use a quadratic polynomial fit, using four 
photometric observations around the spectroscopic date. 

In this study, we restricted the EPM analysis to the first ~ 45-50 days after explosion because 



there a re clear departures from linearity in the 9/v versus t plots after this date. iDessart &: Hilliei 



(|2005bl ) show that the D05 models are poor at late time and therefore should not be used for such 
epochs, supporting this restriction. In Figure[TT]we plot the EPM solutions for SN 1999em (because 
it has extensive photometric and spectroscopic observations during the plateau phase) using the 
{BV}, {BVI}, and {VI} filter subsets and the D05 models. The solid line corresponds to the 
least-squares fit to the derived EPM quantities using the first ~ 70 days after explosion, while the 
dashed line correspond to the least-squares fit using only the first ~ 40 days after explosion. As can 
be noted, after ~ 40 days from explosion (marked with a red triangle) there is departure from the 
linear 6/v versus t relation in all three cases. This justifies our second and third selection criteria 
in § 2.3. However, this restriction severely lower s the number of SNe of our sample to which we can 



apply the EPM. Out of the initial 51 SNe of the lHamuy et al.l (120081 ) sample, only 11 objects fulfill 



the requirement of having a plateau behavior and having early-time photometry and spectroscopy 
for the EPM analysis. 
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3.6.1. EPM Analysis of Individual SNe 

We present here the EPM analysis of 12 SNe IIP (11 from our database and one from the ht- 
erature) with early spectroscopic and photometric observations. We carried out the analysis using 
three different filter subsets {{BV}, {BVI}, {VI}), two sets of host-galaxy extinction (OLI, DES), 
and two atmosphere models (E96, DOS), yielding a total of 12 solutions for each SN. The tables that 
summarize the EPM quantities are available in electronic format for all 144 cases. In the remainder 
of this section we restrict the presentation to the 6 solutions that use the DES extinction because 
they give the lowest dispersion in the Hubble diagrams. Figures [T2H231 show these 6 solutions for 
each of the 12 SNe. Below, we provide the EPM distance D and the explosion date to ^-nd their 
uncertainties, using DES and the {VI} filter subset, and we compare the time of explosion to the 
range restricted by pre-SN images of the host galaxies. These results are also summarized in Table 

El 

In order to obtain a more realistic estimation of the uncertainty in the distance and the ex- 
plosion date, we computed 100 Monte Carlo simulations for each SN, in which we varied all the 
parameters involved in the EPM (see Table [6|), and we averaged the 100 distances and explosion 
dates to derive the EPM values of D and to- This produces small differences between the re- 
sults computed from the initial single EPM solution and that obtained from the 100 Monte Carlo 
simulations, but the latter provides a much more realistic estimate of the uncertainties. 



SN 1992ha 



Figure [121 shows Q/v versus time for SN 1992ba using the {BV}, {BVI}, and {VI} filter subsets 
and the E96 and D05. We used 3 epochs (JD 2448896.9-2448922.8) to compute the distance to this 
SN. In order to use the velocities measured on JD 2448896.9 and 2448900.9, we had to extrapolate 
the /-band photometry until JD 2 448896 .9. 



SN 1992ba was discovered bv lEvansI (|l992l ) on JD 2448896.3. iMcNaughtl (Il992l ) reported that 
the SN was not present on a plate taken on JD 2448883.2 with limiting magnitude 19. The EPM 
solution yields to = JD 2448883.9 ± 3.0 using the E96 models and to = JD 2448879.8 ± 5.6 with 
D05. These results agree (within la) with the explosion date constrained by the pre- and post- 
explosion observations. The distances derived to SN 1992ba are D = 16.4 zt 2.5 Mpc and D = 27.2 
lb 6.5 Mpc using the E96 and the D05 dilution factors, respectively. 



SN 1999hr 

Figure [131 shows 6 /v versus time for SN 1999br using the {BV}., {BVI}., and {VI} filter 
subsets and the D05 atmosphere models. We used 5 epochs (JD 2451291.7-2451309.7) to compute 
the distance to this SN. The EPM solution shows some departure from linearity using the {BV} 
and {BVI} filter subsets. SN 1999br exhibits very low expansion velocities, so we were unable to 
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obtain its distance using the E96 models. This is because the photospheric velocity conversion 
factor Ynp/Vphot is not defined at low expansion velocities (see § 3.4 and Figure [8|). The EPM 
solution yields to = JD 2451275.6 it 7.7 using the DOS models. This result compares very well with 
the obs ervations, because SN 1999br was discovered by the Lick Observatory Supernova Search 
(LOSS; iFilippenko et al.l l200ll ) with the Katzman Automatic Imaging Telescope (KAIT) on JD 
2451280.9 (|Kinglll999l ). An im age tak en on JD 2451264.9 showed nothing at the SN position at a 
limiting magnitudes of 18.5 |Lil ll999al ). The EPM distance to SN 1999br is D = 39.5 ± 13.5 Mpc 
using the D05 dilution factors. 



SN 1999em 



SN 1999em is the best-ever observed SN IIP. Many photometric and spectroscopic observations 
were made by different observers during the plateau phase. Figure O shows 0/v versus time for 
the SN 1999em using the {BV}, {BVI}, and {VI} filter subsets and the E96 and D05 models. 
Table [7] summarizes the EPM quantities derived from the {VI} filter subset. We used 25 epochs 
(JD 2451482. 8 -2451 514.8) to derive the dis tance to SN 1 9 99em . Four spectra were taken from 
Hamuy et al. ( 2001 ) and the other 21 from Ejeonard et al.l (2002). In some cases there were two 
spectra taken at the same epoch from both sources; we used them individualy in the EPM solution 
instead of averaging the measured velocities from each spectrum. We removed the first spectrum 
(JD 2451481.8) from the EPM solution because it shows a clear departure from the linear 9 /v versus 
t relation. The EPM solutions using the E96 and D05 models are quite linear and show great detail 
in the evolution oi 9 /v due to the high-quality spectroscopic and photometric coverage. However, 
the E96 solution shows a small departure from linearity in the last two spectroscopic epochs. This 
effect is probably due to the high rise in the Vup/Vphot ratio at low_velocities in the E96 models. 

SN 1999em was discovered on JD 2451480.9 by the LOSS (iLil ll999bl ). An image taken at the 
position of the SN on JD 2451472.0 showed nothing at a limiting magnitude of 19.0. The EPM yields 
to = JD 2451476.3 ± 1.1 and to = JD 2451474.0 ± 2.0 using the E96 and D05 models, respectively. 
These explosions dates are between the pre-discovery and discovery dates. The distances derived 
to SN 1999em are L> = 9.3 ± 0.5 Mpc from E96 and D = 13.9 ± 1.4 Mpc from D05. 



SN 1999gi 

Figure [15] shows 9/v versus time for SN 1999gi using the {BV}, {BVI}, and {VI} filter 
subsets and the E96 and D05 models. We used 5 epochs (JD 2451525.0- 2451556.9) to apply th e 



EPM method. All of the spectra and the photometry were taken from [Leonard et al.l (|2002bl ). 
The first spectrum (JD 2451522.9) was removed from the EPM solutions because it yields an H/3 
velocity of ~ 26000 km s~^, well above the range of the photospheric velocity conversion (see § 3.4 
and Figure [8]). The explosion dates of SN 1999gi obtained using the EPM are to = JD 2451517.0 
± 1.2 using the E96 models and to = JD 2451515.6 it 2.4 with D05. These results agree with the 
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observations because a pre-discovery image taken on JD 2451515.7 (ITrondal et al.lll999l ) showed 
nothing at the SN position (Umiting unfiltered m agnitude of 18.5). SN 1999gi was discovered on 
JD 2451522.3 (jNakano. Sumoto. Sz Kushidall2002l ) on unfiltered CCD frames, so the explosion date 
can be constrained in a range of only 6.6 days. We derive a distance of -D = 11.7 it 0.8 Mpc and 
D = 17 A lb 2.3 Mpc using the E96 and D05 models, respectively. 



SN 2002gw 



Figure [TU] shows 9/v versus time for SN 2002gw using the {BV}, {BVI}, and {VI} filter 
subsets and the E96 and D05 models. The EPM solutions were obtained using 6 epochs (JD 
2452573.1-2452590.7). The EPM yields explosion times of to = JD 2452557.9 ± 2.7 and to = JD 
2452551.7 it 7. 6 (using E96 and D05 dilution factors, respectively). SN 2002gw was discovered on 
JD 2452560.8 (jMonardlbooa ). An image taken on JD 2452529.6 shows nothing at the SN position 
at a limiting magnitud e of 18.5. Also, an unfilt ered CCD image taken on JD 2452559.1 shows the 
SN at magnitude 18.3 (jitagaki fc Nakand l2002l ) . The EPM explosion dates are in agreement with 
the SN explosion date constrained by the observations. The EPM distances are D = 37.4 it 4.9 
Mpc and D = 63.9 it 17.0 Mpc using E96 and D05, respectively. 



SN 200 3 T 

Figure dZ] shows B/v versus time for SN 2003T using the {BV], {BVI] and {VI] filter subsets 
and the E96 and D05 models. The EPM solutions were obtained using 3 epochs (JD 2452667.9- 
2452701.7). The EPM explosion dates are to = JD 2452654.2 ± using E96 models and to = JD 
2452648.9 it 3.4 with D05. In both cases the third epoch used to derive the distance is beyond 
~ 45 days after the EPM to, but it proves neccesary to include it to compute the EPM analysis. 
This SN was discovered by the Lic k and Tenagra Observatories Supernova Search (LOTOSS) on 



JD 2452664.9 (jSchwartz Lill2003l ). An image taken on JD 2452644.9 shows nothing at a limiting 
magnitude of 19.0, in good agreement with the EPM analysis. The EPM distances are D = 87.8 
± 13.5 Mpc using E96 and D = 147.3 ± 35.7 Mpc with D05. 



SN 2003hl 

Figure [TH] shows 9/v versus time for SN 2003bl using the {BV}^ {BVI}^ and {VI} filter subsets 
and D05 models. The EPM solutions were obtained using 4 epochs (JD 2452701.8-2452735.8). As 
with the SN 1999br, we were unable to apply the EPM using E96 because we only had two spectra 
with velocities higher than 4500 km s^^, and so the photospheric velocity correction could not 
be apphed (see ^ 3.4 an d Figure E]). SN 2003bl was discovered by LOTOSS on JD 2452701.0 
(jSwift. Weisz. &: Lill2003l ). A pre-discovery image taken on JD 2452438.8 shows nothing at the SN 
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position at a limiting magnitud of 19.0. The EPM yields to = JD 2452692.6 ± 2.8, consistent with 
the SN discovery date. The EPM distance is D = 92.4 ± 14.2 Mpc. 



SN 2003hn 

Figure [19] shows O/v versus time for SN 2003bn using the {BV}, {BVI}, and {VI} filter 
subsets and the E96 and D05 models. We computed the EPM analysis using 3 epochs (JD 
2452706.6-2452729.7). The EPM yields explosions dates of to = JD 2452693.4 ± 2.7 and to = 
JD 2452687.0 ± 9.0 from E96 and DOS , respectively. SN 2003bn was discovered on JD 2452698.0 



(jWood-Vasev. Aldering NugentI l2003l ) . Two pre-discovery NEAT images shows nothing at the 
SN position on JD 2452691.5 (limiting magnitude of 21.0) and the SN at a magnitude of 20.2 on 
JD 2452692.8, which restricted the explosion date in a range of only 1.3 days. This value for to is 
in agreement within one a with the EPM to derived using E96 and D05. The EPM distances from 
E96 and D05 are D = 50.2 ± 7.0 Mpc and D = 87.2 ± 28.0 Mpc, respectively. 



SN 2003ef 



Figure [20l shows 9/v versus time for SN 2003ef using the {BV}, {BVI}, and {VI} filter subsets 
and the E96 and D05 models. We computed the EPM analysis using 4 epochs (JD 2452780.7- 
2452797.6). The explosion date derived are to = JD 2452759.8 ± 4.7 and to = JD 2452748.4 ± 15.6 
with E96 and D05, respectivel y. SN 2003ef was discovered by the LOTOSS on JD 2452770.8 (mag. 
about 16.3) (jWeisz Lilbood ). consistent with the EPM value of to. A KAIT image taken on JD 
2452720.8 showed nothing at the SN position at a limiting magnitude of 18.5. The EPM distances 
are D = 38.7 ± 6.53 Mpc with E96 and D = 74.4 ± 30.3 Mpc with D05. 



SN 2003hl 



Figure El] shows 9/v versus time for SN 2003hl using the {BV}, {BVI}, and {VI} filter subsets 
and the E96 and D05 models. The EPM solutions were obtained using 3 epochs (JD 2452879.9- 
2452908.7). We estimated the explosion dates on to = JD 2452872.3 ± 1.7 and to = JD 2452865.4 
lb 5.9 using E96 and D05, respectively. SN 2003hl was discovere d on JD 2452872.0 during the 
LOTOSS program at a magnitude of 16.5 ([Moore. Li. k. Bolesll2003l ). A pre-discovery KAIT image 
taken on JD 2452863.0 shows nothing at the SN position at a limiting magnitude of 19.0. This 
image restricts the explosion date in a range of 9 days. The EPM explosion dates are in agreement 
with the observations (within one a). We derived EPM distances of -D = 17.7 it 2.1 Mpc with E96 
and D = 30.3 ± 6.3 Mpc with D05. 
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SN 2003hn 



Figure [22] shows O/v versus time for SN 2003hn using the {BV}, {BVI}, and {VI} filter 
subsets and the E96 and D05 models. The EPM solutions were obtained using 4 epochs (JD 
2452878.2-2452900.9). The EPM explosion dates derived are to = JD 2452859.5 ± 3.8 and to = 
JD 2452853.8 it 9.3 using the E 96 and 005 d ilution factors, respectively. This SN was discovered 
on JD 2452877.2 at 14.1 mag by lEvand (|2003l ). Evans also reported that the SN was not visible at 
15.5 mag on JD 2452856.5. This date agrees with the explosion date derived from E96 and is < Icr 
lower than that derived from D05. The EPM solutions lead to D = 16.9 it 2.2 Mpc and D = 26.3 
±7.1 Mpc using E96 and D05 models, respectively. 



SN 2003iq 

Figure [23] shows 9/v versus time for SN 2003iq using the {BV}, {BVI}, and {VI} filter subsets 
and the E96 and D05 models. The EP M solutions were obtained using 4 epochs (JD 2452928.7- 



2452948.7). This SN was discovered bv [Llaoasseti (|2003l l on JD 2452921.5, while monitoring SN 



2003hl in the same host galaxy. A pre-discovery image taken on 2452918.5 shows nothing at the 
SN position. These reports constrain the explosion date to a range of only three days. The EPM 
method yields to = JD 2452909.6 ± 4.3 using E96 and to = JD 2452905.6 ± 9.5 using D05. In 
both cases the explosion date is far earlier than expected because the SN was not present on JD 
2452918.5. This implies that the EPM solutions to this SN are not satisfactory. We derived EPM 
distances o{ D = 36.0 ± 5.6 Mpc with Em and D = 53.3 ± 17.1 Mpc with D05. 



4. Discussion 
4.1. External Comparison 
4.1.1. Previous EPM distances 
The EPM met hod has already been applied to SN 1999em by other authors, as follows. (1) 



Hamuy et al.l (|200ll ) employed the E96 dilution factors and eight different filter subsets to perform 
the EPM analysis to this SN, using a cross-correlation technique to estimate the photospheric 
velocity and adopting a host-galaxy extinction of Ay = 0.18 mag. They derived a distance of 6.9 
± 0.1, 7.4 ± 0.1, and 7.3 ± 0.1 Mpc from the {BV}, {BVI}, and {VI} filter subsets, respectively. 
These values are in agreement with our estimates of 6.9 it 0.6, 7.5 it 0.6, and 9.3 it 0.5 Mpc (from the 



BV}, {BVI} and {VI} filter subsets, respectively), except in the {VI} case. (2) [Leonard et al 



BV} 
2OO2I : 



(|2002l ) employed the E96 models to derive the distance to SN 1999em. They used four weak 
unblended spectral features (Fe 11 AA 4629, 5276, 5318, and Sc 11 A 4670) as the photospheric velocity 
indicators. They adopted a host galaxy reddening of Ay = 0.31 mag, the same value predicted by 
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DES. They derived a distance of 7.7 ± 0.2, 8.3 ± 0.2 and 8.8 ± 0.3 Mpc from the {BV}, {BVI} 
and jVI} filter s ubsets , respectively. These results are in agreement with our E96 distances. (3) 



Elmhamdi et al.l (|2003l ) determined a color temperature by fitting blackbody functions to their 
observed spectra of SN 1999em. They used the E96 models and adopted a reddening of Ay = 0.31 
to derive a distance of 7.8 ± 0.3 Mpc to t his SN. This resul t is in agreement with our EPM distance 
using E96, except in the {VI} case. (4) IPessart Hillierl ()2006l ) applied the EPM method to SN 
1999em using E96 and DOS. They adopted the SN 1999em DES reddening value of Ay = 0.31 
mag. Using the E96 models, they derived a distance of 8.6 ± 0.8, 9.7 it 1.0, and 11.7 ± 1.5 Mpc 
from the {-BF}, {BVI}, and {VI} filter subsets, respectively, which are somewhat greater than 
our distances. They also used specific D05 models for this SN to apply the EPM. They found a 
distance of 12.4 ± 1.2, 12.4 ± 1.3, and 12.4 ± 1.3 Mpc from the {BV}, {BVI}, and {VI} filter 
subsets, respectively, using 7 epochs and 11.7 ± 1.0, 11.6 ± 1.0 and 11.5 ± 0.9 Mpc from the {BV}, 
{BVI}, and {VI} filter subsets, respectively, using 8 epochs. These results compare very well with 
our distances using D05 of 11.6 ± 1.2, 12.1 ± 0.9, and 13.9 ± 1.4 Mpc from the {BV}, {BVI}, 
and {VI} filter subsets, respectively. 



4.1.2. SEAM Distance 



SEAM is a technique similar to the EPM, but it avo ids the use of dilu tion factors by doing 
synthetic spectral fitting to the observed spectra of the SN. iBaron et al.l (|2004l ) applied this method 
to SN 1999em. They derived a distance of D = 12.5 it 2.3 Mpc, in good agreement with our 
distances derived using the D05 models (11.6 ± 1.2, 12.1 ± 0.9, and 13.9 it 1.4 Mpc from the 
{BV}, {BVI}, and {VI} filter subsets, respectively), but significantly greater than the EPM 
distances derived using E96 (6.9 ± 0.6, 7. 5 ± 0.6 and 9.3 ± 0.5 Mpc from the {BV}, {BVI} 
and {VI} filter subsets, respectively). Also lDessart fc Hillierl (j2006l ) applied this technique to SN 
1999em, deriving a SEAM distance of 12.2 Mpc (11.5 Mpc) using seven (eight) observations, in 
agreement with the iBaron et al.l ()2004 ) results. 



4.1.3. Cepheid Distance 



Leonard et al.l (|2003l ) identified 41 Cepheid variable stars in NGC 1637, the host galaxy of SN 



1999em. They derived a Cepheid distance to NGC 1637 of D = 11.7 ± 1.0 Mpc. As with the 
SEAM results, the Cepheid distance is consistent with our EPM distances from the D05 models 
(11.2 ± 0.2, 12.0 ± 0.2, and 14.0 ± 0.2 Mpc from the {BV}, {BVI}, and {VI} filter subsets, 
respectively). In all cases, the E96 models lead to significantly lower distances (6.9 it 0.6, 7.5 it 
0.6, and 9.3 ± 0.5 Mpc from the {BV}, {BVI}, and {VI} filter subsets, respectively). 
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4.2. Error Analysis 

4.2.1. Effects of Reddening 



While the ISchlegel et alj (|l998l ) IR maps provide a reasonable estimate of the amount of 
Galactic foreground extinction, the determination of host-galaxy extinction is a more challenging 
task. This is a potential problem because the distances derived using EPM depend on the adopted 
host-galaxy extinction. In order to investigate the sensitivity of the distances to dust extinction, we 
performed the EPM analysis of all the SNe in our sample using the {VI} filter subset by varying 
the amount of host-galaxy visual extinction Ay in steps of IS. Ay = 0.1 mag. Figure [24l shows the 
normalized EPM distances as a function of host-galaxy visual extinction Ay relative to the DES 
value (A Ay = 0). As can be seen, the EPM is quite insensitive to the amount of host galaxy 
extinction adopted. On average, the distances change by less than ~ 10 % from /S.Ay = 0.0 to 
ISAy = 0.5 mag and by less than ~ 20% from ISAy = 0.0 to AAy = —0.5 mag. Therefore, even a 
large systematic error of 0.5 mag in Ay will no t tran slate into a large error in the EPM distance. 
This effect is discussed in Eastman fc Kirshnerl (jl989l ). 



4.2.2. Other Sources of Error 

Table El lists all the error sources in EPM and their typical values. In order to investigate which 
source contributes the most to the uncertainty in the EPM distance, we performed the EPM analysis 
of SN 1999gi (whose photometry and spectroscopy coverage is representative of our sample), and 
we changed the error of a single source (listed in Table [6]) leaving all others unchanged. We found 
two main sources of error. In the E96 case, the errors in the photospheric velocity conversion and 
the dilution factors have the largest effect in the distance uncertainty, each one contributing ~ 30 % 
of the total error, while in the D05 case the error in the dilution factors produces ~ 70 % of the 
uncertainty in the distance, far greater than that due to the error in the photospheric velocity 
conversion (~ 10 % of the total error). All of the other errors have a secondary effect in the total 
error. 



4.3. Hubble Diagrams 



Since the discovery of the expansion of the Universe (Illubblelll929l ). the determination of the 
expansion rate, the Hubble constant (Hq), has become one of the most important challenges in as- 
tronomy and cosmology. Using the velocity-distance relation ( Hubble diagram) calibrated with the 
Cepheid period-luminosity relation. iHubble &: HumasonI (193l|) obtained Hq ~ 500 km s^^ Mpc^^. 
During the second half of the 20th century, the Cepheid relation was significantly improved, and new 
Hubble diagrams were obtained, yielding Hubble constants in the range ~ 50-100 km s~^ Mpc~^. 
Today the discrepancy is not over, but there is a convergence into a value of Hq ~ 65-80 km s~^ Mpc~^) 



-17- 



(|Sandage et alJbood : iFreedman et alJboOll : iRiess et alJboOSi '). 

In this work we applied the EPM method to 12 SNe using two sets of dilution factors (E96, 
D05), two extinction determination methods (OLI, DES), and three filter subsets {{BV}, {BVI} 
and {V^/}) to derive their distances. In order to obtain the host-galaxy redshifts relative to the 
cosmic microwave background (CMB), we corrected the heliocentric host-galaxy redshifts for the 
peculiar velocity of the Sun relative to the CMB rest frame. F or this purpose we added a velocity 
vector of 371 km s~^ in the direction {l,b) = (264.14°, 48.26°) (jpixsen et al.lll99d ) to the heliocen- 
tric redshifts. The resulting CMB redshifts are given in Table [1] 

Using the CMB host-galaxy redshifts, we constructed 12 different Hubble diagrams. In each 
case we computed a linear fit weighting the error in distance and redshift (assumed to be 300 km 
for all SNe) in order to derive the Hubble constant. Figures [25VI27I show the Hubble diagrams ob- 
tained with OLI reddenings, from the {BV}, {BVI} and {VI} filter subsets, respectively. Figures 
[28H30] show the same diagrams but this time using DES extinctions. Each diagram is labeled with 
the derived Hubble constant, the reduced and the dispersion in distance modulus from the 
linear fit. The resulting Hq values are summarized in Table El 

There is a systematic difference in the Hq values obtained with the E96 and DOS models. Using 
E96 we obtained Hq = 89-101 km s~^ Mpc~^ while D05 yielded Hq = 52-66 km s"^ Mpc"^ This 
difference arises both from the systematically higher D05 dilution factors and the different photo- 
spheric velocity conversion between both models. These two effects combined lead to differences of 
~ 40 % in the EPM distances. 

The use of different filter subsets leads to Hq values consistent within 1 a for a fixed atmosphere 
model. This is a very important result, because it suggest an internal consistency for each set of 
atmosphere models. However, the use of different filter subsets produces significant differences in 
dispersion, increasing from ~ 0.3 ({VI}) to ~ 0.4 ({BVI}) and o"^ 0.5 ({BV}) (see Table 



The special case of D05 with {VI} and DES, leads to = 0.32, which corresponds to ^ 15 



70 



of error in distance. Clearly when the B band is employed, the dispersion in the Hubble diagram 
increases considerably. This is due to the fact that the dilution factors in the B band, are mostly 
determined by the effect of line blanketing rather than by electron scattering, so the assumption 
of a "dilute" Black Body becomes less reliable. It is also possible that metallicity differences of 
the SNe could explain part of this scatter, although both atmosphere model sets predict a modest 
effect of metallicity in the emergent flux at wavelengths longer than ~4000 A. 

As expected, it can be noted that there are no significant differences in the Hq values and in 
the Hubble diagram dispersion between the DES and OLI reddening methods. This is because 
there is no systematic difference in the reddening between both methods (see § 3.5). However, the 
DES method leads to somewhat lower dispersion in the Hubble diagrams than the OLI technique. 

Finally, SN 2003hl and SN 2003iq are of particular interest because they both exploded in the 
same galaxy. To our disappointment, all 12 possible combinations of filter subsets, reddening, and 
atmosphere models lead to significant differences in the EPM distance to the host galaxy. The most 
extreme case is the {BV}, E96, and OLI combination, which leads to a distance of 32.5 it 8.5 Mpc 
to SN 2003iq and 12.8 ± 1.6 Mpc to SN 2003hl (a difference of 2.3 sigma). The smallest discrepancy 
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occurs with the {VI}, D05, and DES combinations (30.3 ± 6.3 and 53.3 ± 17.1 Mpc for SN 2003hl 
and SN 2003iq, respectively), which is also the combination that produces the lowest dispersion in 
the Hubble diagram. As discussed in §3, the EPM solutions to SN 2003iq yield an explosion time 
inconsistent with a pre-discovery image, hence the EPM distance to SN 2003iq is quite suspicious. 
The light curve and the spectroscopic velocity evolution of both SNe are typical of the Type IIP 
class, so this is not the reason for the inconsistency in the EPM solutions. Given the high host- 
galaxy extinction to SN 2003hl, it is conceivable that the discrepancy could be due to a departure 
from the standard reddening law. We explored this possibility by varying Ry over the range 1.1- 
3.1 using the {VI} and D05 models. For the case i?v'=3.1 we obtain D(03hl)/D(03iq)=0.57, while 
for Rv=l-1 we get D(03hl)/D(03iq)=0.72. This exercise shows that, although the discrepancy is 
not completely remo ved, lowering Ry brings the two distances in better agreement. Interestingly, 



Olivares et al.l (j2008l ) found an optimal value of iiy=1.4 from a set of 34 SNe IIP using the SCM. 



All this suggests that either 1) the dust around SNe IIP has different properties than the Galactic 
dust, or 2) SNe IIP are surrounde d by Galactic-l ike dust whose geometric distribution is responsible 
for an abnormal low value of Ry (|Goobaiil2008l ). 



4.4. External Calibration and the Internal precision of the EPM 

In the previous section we have shown that there is a systematic difference in the Hq values 
derived using the E96 and the DOS models. In order to remove this systematic effect, we applied a 
calibration factor (given by the ratio between some external Hq value and the EPM Hq value) to the 
distances derived using E96 and D05. For this purpose we u sed the value of Hq = 72 km s^^ Mpc~^ 



derived from the HST Key Project (jFreedman et al.ll200ll ). This external calibration allows us to 



bring the EPM distances to the Cepheid scale and allows us to remove the systematic difference in 
the EPM distances between E96 and D05. Figure [3T] shows (top panel) the DOS distances versus 
the E96 distances divided by a calibration factor of 1.37 and 0.79, respectively. In both cases the 
EPM distances were derived using the {VI} filter subset and the DES reddening. As can be seen, 
after applying this correction, the systematic differences disappear. The dashed line in the top 
panel corresponds to the one to one relation. Also, in Figure [31] (bottom panel) are plotted the 
differences between the corrected E96 and DOS distances, normalized to the corresponding average 
between the corrected E96 and DOS distance. We found a standard deviation of o" = 0.12. Since the 
dispersion arises from the combined errors in the E96 and the DOS distances, the internal random 
errors in any of the EPM implementation must be less than 12%. Note that this scatter is smaller 
than the ~ 15 % dispersion seen in the Hubble diagrams, which is affected by the peculiar motion 
of the host galaxies. The 12% scatter is independent of the redshift and must be an upper value of 
the internal precision of the EPM. 
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5. Conclusions 

In this work we iiave applied the EPM method to 12 SNe IIP. We contructed 12 different Hub- 
ble diagrams, using three different filter subsets {{BV}, {BVI}, {VI}), two atmosphere models 
(E96, DOS), and two methods to determine the amount of host-galaxy extinction (DES, OLI). Our 
main conclusions are as follows. 

The EPM must be restricted to the first ~ 45-50 days from explosion. After that epoch the 
method may display departures from linearity in the 9/v versus time relation, and therefore an 
internal inconsistency. 

The results are less precise when the B band is used in the EPM analysis, regardless of the 
atmosphere models employed (E96 or D05). The dispersion in the Hubble diagrams increases con- 
siderably from 0.3 to 0.5 mag when the B band is included and the V band is removed from the 
filter subset. Despite the loss in precision, there are no significant differences in the resulting dis- 
tances. 

We investigated the effect of host-galaxy reddening in the EPM distances. For this purpose 
we computed many EPM solutions varying the amount of visual extinction, and we found that a 
difference of AAy = 0.5 mag leads on average to a difference of ~ 5-10 % in distance. Therefore, 
we conclude that the method is quite insensitive to the effect of dust. 

Systematic differences in the atmosphere models lead to ~ 50 % differences in the EPM dis- 
tances and to values of Hq between 52 and 101 km Mpc~^. This effect is due to the systematic 
difference in the photospheric velocity conversion and the dilution factors. The latter is currently 
the greatest source of uncertainty in the EPM method. 

The Hubble diagram with the lowest dispersion (cr^ = 0.32 mag) was obtained using the com- 
bination D05, {^/}, DES. Despite the systematic uncertainties in the EPM, this dispersion is quite 
low and corresponds to a prec ision of ^ 15 % in dis t ance. This precision is similar to that of the 



SCM method for type II SNe (jHamuy k Pintoll2002l : lOlivares et al.ll2008l) an d to the Tully-Fisher 



relation for spiral galaxies with a dispersion of o" ~ 0.30 mag (ISakai et al.ll2000l ). However, the EPM 
dispersion is greater than that of the M/Ami5 relation for Type la SNe, which has a dispersion of 
a ~ 0.15-0.20 mag; however, if the EPM is applied to a sample of SNe IIP in the Hubble flow, the 
dispersion in the Hubble diagram might decrease. 

Finally, despite the systematic differences in the Hq value, EPM has great potential as an 
extragalactic distance indicator and can potentially be applied to a sample of high-redshift SNe IIP 
in order to check in an independent way the accelerating expansion of the universe. 
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Fig. 1. — Optical light curves of four SNe during the first ~ 120 days of their evolution. The top 
axis of each panel gives the phase in days since the EPM explosion time, derived using the DOS 
models (see Table [5]) . 
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Fig. 2. — Optical light curves of four SNe during the first ~ 120 days of their evolution. The top 
axis of each panel gives the phase in days since the EPM explosion time, derived using the DOS 
models (see Table [5]) . 
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Fig. 3. — Optical light curves of four SNe during the first ~ 120 days of their evolution. The top 
axis of each panel gives the phase in days since the EPM explosion time, derived using the DOS 
models (see Table [5]) . 
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Fig. 4. — Dilution factors (" as a function of the color temperature, computed at z = from 
the E96 (blue dots) and D05 (red dots) atmosphere models for three different filter subsets 
{{BV}, {BVI}, {VI}). The blue (red) line corresponds to the polynomial fit performed to the 
E96 (DOS) dilution factors. 
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Fig. 5. — Line velocity evolution determined from the P Cygni absorption minima of four different 
features during ~ 100 days after discovery. The top axis of each panel gives the phase in days since 
the EPM explosion time derived using the DOS models (see Table [5]) . 
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Fig. 6. — Line velocity evolution determined from the P Cygni absorption minima of four different 
features during ~ 100 days after discovery. The top axis of each panel gives the phase in days since 
the EPM explosion time derived using the DOS models (see Table [5]) . 
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Fig. 7. — Line velocity evolution determined from the P Cygni absorption minima of four different 
features during ~ 100 days after discovery. The top axis of each panel gives the phase in days since 
the EPM explosion time derived using the DOS models (see Table [5]) . 



- 30 - 




Fig. 8. — Ratio of the H/3 velocity to the photospheric velocity versus the H/J velocity of the 
individual SN models. The blue dots correspond to E96 models and the red dots to D05 models. 
The blue (red) line corresponds to the polynomial fit performed to the E96 (D05) photospheric 
velocity conversion. 
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Fig. 9. — Ratio of the Ha velocity to the H/3 velocity as a function of the H/3 velocity. The triangles 
and the squares represent velocities measured from the spectra of our SN sample. The open and 
filled black circles correspond to the velocity ratio measured from the synthetic spectra of E96 and 
D05, respectively. 
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Fig. 10. — Comparison between the DES and OLI reddening methods for the 12 SNe. The straight 
hne has a slope of unity. The more deviant SNe are labeled. 
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Fig. 11. — The ratio 9/v as a function of time for SN 1999em using the {BV}, {BVI}, and {VI} 
filter subsets and the D05 models. The solid (dashed) lines correspond to unweighted least-squares 
fits to the derived EPM quantities using ~ 70 (40) days after explosion. The red triangle in the 
bottom panel shows day ~ 40 after explosion. 



- 34 - 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




880 890 900 910 920 

JD-2448000 



Fig. 12.— The ratio Ojv as a function of time for SN 1992ba using the {BY}, {BVI], and {VI} 
filter subsets. The ridge lines correspond to unweighted least-squares fits to the derived EPM 
quantities. The upper and lower panels show the results using E96 and DOS dilution factors, 
respectively. In all cases we employ the DES reddening. 



- 35 - 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




270 280 290 300 310 



JD-2451000 



Fig. 13.— The ratio O/v as a function of time for SN 1999br using the {BV], {BVI}, and {VI] 
fiher subsets and the DOS models. The ridge hues correspond to unweighted least-squares fits to 
the derived EPM quantities. In all cases we employ the DES reddening. 
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Fig. 14. — The ratio 9/v as a function of time for SN 1999em using the {BV}^ {BVI}, and 
{VI} filter subsets. The ridge hnes correspond to unweighted least-squares fits to the derived 
EPM quantities. The upper and lower panels show the results using E96 and DOS dilution factors, 
respectively. In all cases we employ the DES reddening. 
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Fig. 15.— The ratio Qjv as a function of time for SN 1999gi using the {BV}, {BVl}, and {VI} filter 
subsets. The ridge hnes correspond to unweighted least-squares fits to the derived EPM quantities. 
The upper and lower panels show the results using E96 and D05 dilution factors, respectively. In 
all cases we employ the DES reddening. 
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Fig. 16. — The ratio 9/v as a function of time for SN 2002gw using the {BV}, {BVI}, and 
{VI} filter subsets. The ridge hues correspond to unweighted least-squares fits to the derived 
EPM quantities. The upper and lower panels show the results using E96 and DOS dilution factors, 
respectively. In all cases we employ the DES reddening. 
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Fig. 17.— The ratio Ojv as a function of time for SN 2003T using the {BV}, {BVl], and {VI} filter 
subsets. The ridge fines correspond to unweighted least-squares fits to the derived EPM quantities. 
The upper and lower panels show the results using E96 and D05 dilution factors, respectively. In 
all cases we employ the DES reddening. 
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Fig. 18.— The ratio O/v as a function of time for SN 2003bl using the {BV}, {BVI}, and {VI} 
filter subsets and the DOS models. The ridge lines correspond to unweighted least-squares fits to 
the derived EPM quantities. In all cases we employ the DES reddening. 
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Fig. 19. — The ratio 0/v as a function of time for SN 2003bn using the {BV}^ {BVI}^ and 
{VI} filter subsets. The ridge hues correspond to unweighted least-squares fits to the derived 
EPM quantities. The upper and lower panels show the results using E96 and DOS dilution factors, 
respectively. In all cases we employ the DES reddening. 
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Fig. 20.— The ratio e/vasa function of time for SN 2003ef using the {BV}, {BVI}, and {VI} filter 
subsets. The ridge hnes correspond to unweighted least-squares fits to the derived EPM quantities. 
The upper and lower panels show the results using E96 and D05 dilution factors, respectively. In 
all cases we employ the DES reddening. 
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Fig. 21.— The ratio O/vasa function of time for SN 2003hl using the {BV}, {BVI}, and {VI} filter 
subsets. The ridge hues correspond to unweighted least-squares fits to the derived EPM quantities. 
The upper and lower panels show the results using E96 and D05 dilution factors, respectively. In 
all cases we employ the DES reddening. 
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Fig. 22.— The ratio Qjv as a function of time for SN 2003hn using the {BV}, {BVl}, and 
{y/} filter subsets. The ridge hues correspond to unweighted least-squares fits to the derived 
EPM quantities. The upper and lower panels show the results using E96 and DOS dilution factors, 
respectively. In all cases we employ the DES reddening. 



-45 - 



— I — < — , — I — < — I — I — I — I — < — \ — I — < — I — I — ] — I — I — I — I — I — ^ 

SN 2003iq 




910 920 930 940 950 

JD-2452000 



Fig. 23.— The ratio O/vasa function of time for SN 2003iq using the {BV}, {BVI}, and {VI} filter 
subsets. The ridge fines correspond to unweighted least-squares fits to the derived EPM quantities. 
The upper and lower panels show the results using E96 and D05 dilution factors, respectively. In 
all cases we employ the DES reddening. 
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Fig. 24. — Normalized EPM distances as a function of the host-galaxy visual extinction relative to 
the DES value (A^y = 0). 
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Fig. 25. — Hubble diagram using the {BV} filter subset and OLI reddening. 
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Fig. 26. — Hubble diagram using the {BVI} filter subset and OLI reddening. 
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Fig. 27. — Hubble diagram using the {VI} filter subset and OLI reddening. 
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Fig. 28. — Hubble diagram using the {BV} filter subset and DES reddening. 
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Fig. 29. — Hubble diagram using the {BVI} filter subset and DES reddening. 
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Fig. 30. — Hubble diagram using the {VI} filter subset and DES reddening. 
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Fig. 31. — Top panel: DOS distances versus E96 distances corrected to the EST Key Project 
Cepheid scale. The dashed line shows slope unity. Bottom panel: differences between the corrected 
distances normalized to the average of the E96 and D05 corrected distances. The 12% scatter 
reflects the internal precision of the EPM. 
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Table 1. Heliocentric and CMB redshifts for the SNe used in this work. 



SN 


Host Galaxy 




source'^ 


CZCMB 






(km s^^) 




(km s^"*^) 


1992ba 


NGC 2082 


1092 


here 


1245 


1999br 


NGC 4900 


960 


NED 


1285 


1999em 


NGC 1637 


800 


L02 


670 


1999gi 


NGC 3184 


543 


here 


831 


2002gw 


NGC 0922 


3117 


here 


2877 


2003T 


UGC 04864 


8368 


NED 


8662 


2003bl 


NGC 5374 


4382 


NED 


4652 


2003bn 


2MASX J10023529-2110531 


3829 


NED 


4173 


2003ef 


NGC 4708 


4440 


here 


4503 


2003hl 


NGC 0772 


2265 


here 


2198 


2003hn 


NGC 1448 


1347 


here 


1102 


2003iq 


NGC 0772 


2364 


here 


2198 



^ The NED values correspond to the redshifts of the host-galaxy nucleus, while 
the values measured in this work ("here") were measured from narrow emission 
lines of H 11 regions at the SN position. Also, L02 corresponds to the value adopted 
from [Leonard et aL ( 2002 ). 
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Table 2. Dilution factor coefficients and dispersion 







E96 






D05 




Filter subset 


bo 


bi 


b2 


a 


bo 


bi 


b2 


a 


{BV} 


0.756 


-0.900 


0.520 


0.048 


0.593 


-0.450 


0.403 


0.075 


{BVI} 


0.733 


-0.693 


0.373 


0.027 


0.711 


-0.476 


0.308 


0.068 


{VI} 


0.702 


-0.531 


0.265 


0.029 


0.915 


-0.747 


0.371 


0.077 



Table 3. H/3 to photospheric velocity ratio coefficients and dispersion. 



j aj(E96) aj(D05) 



1.775 1.014 

1 -1.435 X 10-^ 4.764 x IQ-^ 

2 6.523 X lO-'^ -7.015 x IQ-^ 
cr 0.06 0.04 
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Table 4. SN host galaxy and Galactic extinction adopted. 



SN 


Av (OLI) Av (DES) 


^ Av (IR maps) " 




Host Host 


Galactic 



1992ba 


0.30 


(0.15) 





43 


(0 


16) 





193 


(0 


031) 


1999br 


0.94 


(0.20) 





25 


(0 


16) 





078 


(0 


012) 


1999em 


0.24 


(0.14) 





31 


(0 


16) 





134 


(0 


021) 


1999gi 


1.02 


(0.15) 





56 


(0 


16) 





055 


(0 


009) 


2002gw 


0.18 


(0.16) 





40 


(0 


19) 





065 


(0 


010) 


2003T 


0.35 


(0.15) 





53 


(0 


31) 





104 


(0 


017) 


2003bl 


0.26 


(0.15) 





00 


(0 


16) 





090 


(0 


014) 


2003bn 


-0.04 


(0.15) 





09 


(0 


16) 





215 


(0 


034) 


2003ef 


0.98 


(0.15) 


1 


24 


(0 


25) 





153 


(0 


024) 


2003hl 


1.72 


(0.18) 


1 


24 


(0 


25) 





241 


(0 


039) 


2003hn 


0.46 


(0.14) 





59 


(0 


25) 





047 


(0 


008) 


2003iq 


0.25 


(0.16) 





37 


(0 


16) 





241 


(0 


039) 



Olivares et Zl |2008h . 



Dessart fc Hillieij |2006h ; IPessart et al.l |2008h . 



Schlegel et all (jlQQSh 
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Tablo 5. EPM distances. 



SN 


^(£96) 


io{E96) 




*o(i:»05) 




(Mpc) 


(JD-2448000) 


(Mpc) 


(JD-2448000) 


1992ba 


16.4 (2.5) 


883.1 


9 (3.0) 


27.2 (6.5) 


879.8 (5.6) 


1999br 








39.5 (13.5) 


3275.6 (7.7) 


1999em 


9.3 (0.5) 


3476, 


,3 (1.1) 


13.9 (1.4) 


3474.0 (2.0) 


1999gi 


11.7 (0.8) 


3517, 


.0 (1.2) 


17.4 (2.3) 


3515.6 (2.4) 


2002gw 


37.4 (4.9) 


4557, 


.9 (2.7) 


63.9 (17.0) 


4551.7 (7.6) 


2003T 


87.8 (13.5) 


4654, 


.2 (2.7) 


147.3 (35.7) 


4648.9 (6.1) 


2003bl 








92.4 (14.2) 


4694.5 (2.0) 


2003bn 


50.2 (7.0) 


4693, 


,4 (2.7) 


87.2 (28.0) 


4687.0 (9.0) 


2003ef 


38.7 (6.5) 


4759, 


.8 (4.7) 


74.4 (30.3) 


4748.4 (15.6) 


2003hl 


17.7 (2.1) 


4872, 


.3 (1.7) 


30.3 (6.3) 


4865.4 (5.9) 


2003hn 


16.9 (2.2) 


4859, 


.5 (3.8) 


26.3 (7.1) 


4853.8 (9.3) 


2003iq 


36.0 (5.6) 


4909, 


.6 (4.3) 


53.3 (17.1) 


4905.6 (9.5) 



Note. — The distances were derived using the {VI} filter subset and 
DES reddening. 



Table 6. Error Sources 



Error Source 



Typical Error 



Photometry 
SN redshift 
Foreground extinction 
Host galaxy extinction 
Line expansion velocity 
Photospheric velocity conversion 
Dilution Factors 



0.02 mag 
50 / 200 (km s"^) 
0.02 mag 
0.15 mag 
85 (km s"^) 
0.06 / 0.04 b 
0.03 / 0.07 ^ 



^ Corresponds to the redshifts measured in this work and 
those taken from NED, respectively. 

^ Corresponds to the E96 and DOS models, respectively. 
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Table 7. SN 1999em EPM quantities. 



JD- 


Tvi 


eCvi 


Cvi 


Vphot 


e/vel 


2451000 


(K) 


\ l_U ( llJ 






(km s-i) 


(100 s M|x--^) 


482.8 




( ARQ\ 


0.0321 


(0.0010) 


0.574 


11022 


(JUL). / 


yi 0.\J } 


483.8 






0.0331 


(0.0011) 


0.572 


10355 


ooy.D 


/'81 

( Oi.U ) 


484.8 




yooZ, ) 


0.0341 


(0.0009) 


0.568 


9867 


DUo.O 




485.2 


13810 


(415) 


0.0345 


(0.0011) 


0.566 


9117 




(Q7 7\ 


485.7 


13550 


(414) 


0.0352 


(0.0011) 


0.563 


8942 




f1 09 Pi\ 


485.7 


13544 


(414) 


0.0352 


(0.0011) 


0.563 


8915 


709 1 


luo.u j 


485.8 


13479 


(456) 


0.0353 


(0.0012) 


0.562 


9311 




I yc/. / j 


486.8 


12812 


(425) 


0.0373 


(0.0013) 


0.555 


8817 


t uz.u 




487.9 


11985 


(333) 


0.0403 


(0.0013) 


0.547 


8584 


oO t .0 


(^ 98 (x\ 


488.8 


11587 


(310) 


0.0413 


(0.0013) 


0.544 


8598 


889 8 
OoZ.o 


(^ qq qN 


489.8 


11352 


(256) 


0.0424 


(0.0011) 


0.542 


8476 


921.1 


(138.7) 


491.1 


11077 


(350) 


0.0443 


(0.0016) 


0.541 


7870 


1040.1 


(159.5) 


491.2 


11055 


(358) 


0.0444 


(0.0017) 


0.541 


7824 


1050.6 


(161.4) 


491.7 


10939 


(372) 


0.0453 


(0.0018) 


0.540 


7964 


1053.3 


(162.6) 


492.1 


10840 


(349) 


0.0460 


(0.0018) 


0.539 


7863 


1083.7 


(166.9) 


496.2 


10264 


(312) 


0.0495 


(0.0019) 


0.537 


7031 


1311.6 


(202.6) 


496.7 


10224 


(301) 


0.0497 


(0.0018) 


0.537 


7172 


1290.7 


(199.0) 


501.2 


9610 (224) 


0.0526 


(0.0016) 


0.537 


5921 


1653.2 


(252.9) 


501.7 


9386 (185) 


0.0548 


(0.0014) 


0.538 


6107 


1667.6 


(253.4) 


501.7 


9384 (185) 


0.0548 


(0.0014) 


0.538 


6250 


1630.3 


(247.7) 


501.8 


9362 (189) 


0.0551 


(0.0015) 


0.538 


6506 


1572.5 


(238.9) 


504.8 


8907 (173) 


0.0589 


(0.0016) 


0.542 


5991 


1813.0 


(274.0) 


506.8 


8655 (162) 


0.0605 


(0.0016) 


0.545 


5691 


1950.5 


(293.2) 


510.8 


8248 (63) 


0.0649 


(0.0007) 


0.553 


5156 


2276.6 


(334.0) 


514.8 


7819 (92) 


0.0705 


(0.0013) 


0.565 


4877 


2557.8 


(370.3) 



Note. — The EPM quantities were derived using the {VI} filter subset, the DES 
reddening and the D05 models. 
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Table 8. Summary of Hq values. 





{BV} 


{BVI} 


{VI} 


E96/OLI 


98.4 (9.2) 


100.8 (8.0) 


89.1 (6.9) 


E96/DES 


97.2 (8.7) 


100.5 (8.4) 


91.2 (6.7) 


D05/OLI 


66.2 (4.2) 


60.4 (4.1) 


53.9 (4.3) 


D05/DES 


63.8 (3.9) 


59.6 (4.2) 


52.4 (4.3) 
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Table 9. Summary of dispersions in Hubble diagrams. 





{BV} 


{BVI} 


{VI} 


E96/OLI 


0.53 


0.43 


0.34 


E96/DES 


0.50 


0.41 


0.37 


D05/OLI 


0.57 


0.39 


0.36 


D05/DES 


0.51 


0.37 


0.32 



